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Abstract 

Tectonic patterns on Eutupa are influenced by tidal stress. An important well 
recognized ooniponont. is associated with the orbital erccntricily, which produces a 
diurnaiiy varying stress as Jupiter's apparent position in Europe's sky oscillates in 
longitude. Cycloidal lineaments seem to have formed as cracks propagated in this 
ditirnaliy varying stress field. Maps of theoretical cycloid patterns capture many of 
the characteristics of the observed distribution on Europa. However, a few details 
of the observed cycloid distribution have not been reproduced by previous models. 
Recently, it has been shown that Europa Irks a finite forced obliquity, so Jupiter’s 
apparent position in Europa sky will also oscillate in latitude. We explore this now 
type of diurnal effect on cycloid formation. We find that stress from obliquity may be 
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the key to explaining several dial act eristics of observed cycloids such as the shape 
of eqnaf or- crossing cycloids and the shift in the crack patterns in the Argadmd 
Regies region. Aii of these improvements of the lit between olnervaffoit and theory 
seem to require Jupiter crossing Europe's equatorial plum; 45 c ' to 180'’ after perijove 
passage, suggestive of complex orbital dynamics that Jock the direction of Europe's 
pcricenter with the direction of the ascending node at the Lime these cracks were 
formed. 
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I Introduction 

After the Voyager spacecraft revealed cycloids in the limited area it imaged, 
the Formation of cycloids remained unexplained until Hoppa and Tufl.s (1999) 
proposed that cycloidal cracks form as a result of tensile cracking m response 
to diunmliy varying tidal stresses produced by Europa’s orbital eccentricity. 
In their model, cycloid-shaped cracks form in the following way: When the 
tensile strength of ice is readied, a crack forms perpendicular to the local 
direction of the tensile stress and begins to propagate across the surface. Be 
cause diurnal tidal stress changes in both magnit ude and orientation, cracks 
propagate across an ever-changing stress field. Thus, propagation can follow 
a curving path until it reaches a place and time where the tensile stress is 
insufficient to continue the propagation. The propagation remains dormant 
until later m the orbit, when the stress at the end of the crack once again 
exceeds the ice strength. At that time, propagation continues in a direction 
perpendicular to the new orientation, of tensile stress. The shape of the crack 
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2 The Globa] Picture of Cycloid Formation 

II I Comparison of theoretical models with observation: zero obliqmt/ij 

Fig. 1 shows maps of idealized cycloids that would form at locations across 
Europe's surface, assuming only diurnal stress variations caused by the eccen- 
tricity. The crack shapes shown are those that, propagate from initial erode 
sites at various locations, selected ret intervals of 15° in longitude and at HP 
intervals in latitude between 10° and 70° in the north and south. We also 
model cycloids that form at 5“ latitude both north and south of the equator 
Fig. la. shows crack patterns assuming: a net westward component of prop- 
agation. while Fig. lb assumes the net. sense of propagation is towards the 
east. The tracking patterns depend on the rate of crack propagation and the 
assumed ke strength values (the tension required to start and maintain crack 
propagation). Here we use parameters .similar to those inferred from fitting 
cycloids in the southern hemisphere (Hurford et at, 2007). Specifically, the 
initiation stress is 68 kPa, the propagation strength is 38 kPa arid the effec- 
tive propagation speed is 4.8 km /hr. The crack maps in Fig. i are idealised, of 
course, because details of the local parameters controlling crack propagation 
may influence a specific cycloid 's development . 

For comparison, Fig. 2 shows mapped cycloids actually observed on Europe. 
Bart et al (2003) and Greenberg et at. (2003) noted that hypothetical mod 
tiled cycloids reproduce many of the characteristics of observed cycloids: 1 1 ) 
The long east -west trending cycloids observed in the southern, mid-latitudes 
near longitude .180° (Fig, 2) are similar to hypothetical cycloids at the same 
latitudes, albeit not. the same longitudes (Fig. la&b). In fact, these include the 


cycloids modeled by Hoppa ct al. (2001) and Hurford et al. (2007). (2) Long 
east-west trending cycloid chains run just north and south of the equator in 
Fig. I, as observed near 0° and 180' iu Fig. 2. Such cycloids are prominent 
near 0“ and 180° in Fig. 2 and are not seen at other longitudes because image 
resolution or coverage is inadequate or cycloid features have been disrupted 
by chaotic terrain. (3) Near 90“ longitude in the fat north (Fig, 2) are cycloids 
with north-south orientations. Like the ones in the hypothetical map f Fig. | l. 
they tend to be relatively short and are located far from the equator. (4) The 
formation of angular or “boxy” crack patterns observed near the equator (Fig. 
2} is similar to the model cracks in Fig. 1 in this region. These observed “boxy” 
cycloids lie in Argadnel Regie (often called the “Wedges" region), which is lo 
r ated just south of the equator and west of 180° longitude. The cracks in this 
region, besides demonstrating the boxy cycloid shapes, are often dilated into 
wedge shapes. Similar crack patterns diametrically opposite (in the northern 
hemisphere near IF in longitude) have been noted (Greenberg cl al., 2003; 
Greenberg. 20(15), but there they have not been modified by dilation. 


Besides matching the general characteristics of observed cycloids, maps of 
idealised cycloids might match patterns observed in evdoid distributions. (A) 
In F ig. 1 small zones in which no cycloids form are predicted near (F and 180° 
longitude. In these zones the maximum tensile stress never reaches a level 
high enough to initiate crack formation, thus no cycloids originate from those 
regions. (B) Cycloid patterns in the northern and southern hemisphere arc 
symmetric. In tact, the pattern in each hemisphere is a mirror image of the 
other, reflected about the equator. 
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3 Obliquity arid Cycloid Modeling 


cling cycloids with the inclusion of non synchronous rotation still produces 
zones in which no cycloids form near 0" and 180° longitude but the size of 
these zone increase and their centra's shift: eastward from 180“- However, with 
the reorientation of the surface through iron-synchronous rotation, these zones 
where cycloids fail to form may not be preserved in the observed cycloid die 
tributkm. Most importantly, the distribution of idealized cycloids maintains 
the symmetry between the northern and southern hemisphere. Since idealized 
cycloids that formed with the stress from non -synchronous rotation included 
reproduce the overall characteristics of the observed cycloids, which were also 
reproduced by the diurnal stresses alone and fail to reconcile the detailed mis- 
match between the idealized and observed distribution, we will exclude the 
effect; of stress from non- synchronous rotation on cycloid growth patterns in 
this paper. 

Similarly, a small amount of reorientation hv polar wander may account for 
the mismatch between the model and the observations of cracks. Sat id et nl. 
(2002) found a mismatch between the model of strike-slip displacement and 
observations strike-slip displacement, which they attributed to polar wander 
of Europa’s ice shell. However, polar wander of the type they describe would 
have shifted the Argadnel Regie region north after its formation, implying that 
it formed even further south of its current location, enhancing the asymmetry 
across the equator. 

In summary, the main disagreement is that the observations do not St the 
predicted symmetry across the equator'. Cycloids that cross the equator are 
trot symmetrical, and the "boxy-' cracks just west of 180° are south of the 
equator, while those just west of 0“ are north of the equator. 


The key to modeling the observed characteristics that do not fit the previous 
model (Tig. 1) may be the satellite’s obliquity, c. If Europa’s orbital plane 
were fixed in orientation, then tidal damping would have rapidly removed any 
initial obliquity (Peak, 1974: Gladntan et <U., 1996). However, the orbit plane 
is not fixed, but processes about an axis very close to Jupiter’s spin axis. This 
precession is mainly caused by torques from Jupiter’s oblate figure, but also 
has contributions from interactions with In, Ganvmedo, and Callisto (Lieske, 
.1998). 

As a result of the orbital precession, the configuration to which the spin pule 
is driven by tidal damping is not one of zero obliquity. Instead, as for the 
Moon, the obliquity adjusts to a value such that the spin pole, orbit pule, and 
local invariable pole all remain copianar. Such configurations have beeu culled 
generalized Cassini states (Colombo, 1966; Peaks. 1969: Ward, 1975), because 
they are similar to behavior of the Moon, as first enunciated by G.I). Cassini 
in 1693. 

ff the orbit precession is steady, as assumed by Bills and Ray (2000) in their 
analysis of obliquity for the Galilean satellites, then obliquities would bo tidally 
damped to constant, non-zero values, if the orbit-plane precession rate is pe- 
riodically modulated bv interactions among the neighboring satellites, then 
the tidally damped obliquities will correspondingly fluctuate (Ward and de 
Campli. 1979: Bills, 2005). Given that effect, Europa’s mean obliquity is close 
to O.T. and oscillates by roughly 1(1%. with dominant periods of decades to 
centuries. On still longer time scales (> if! 1 years) the eccentricity and inch- 
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changes in a and 0 while the obliquity causes changes only in It. Jupiter s angu- 
lar position with respect to a fixed location above Europe's surface oscillates 
with an amplitude of 2c: east and west and c: north and south, changing slightly 
the angular distance to the tidal bulge 0. The change in 0 pins daily changes 
in a comhine, yielding the dhirnaily varying part of the tide, producing elastic 
stress on the surface {Sett Appendix A). 


For a water-fax crust, following Hoppa el al. (1989b, 2001) we use the con- 
ventional values of the elastic parameters ft -- 3-52xl() fl Pa and v 0,33 
{Gammon et at.. 1983), and assume a tidal response given by h t — 1.275. 
which corresponds to a diurnal tidal amplitude of changes in the radius of 
30 in. The. coefficient 'AM /Sup,,,, a 3 is approximately 2.5x10 h 


To simulate the formation of a cycloidal crack at a given location, our model 
also requires values for: the effective rate of crack propagation, which largely 
governs length and curvature of the arcs; the sense {eastward or westward) of 
the propagation, which determines the direction of curvature of each arc: the 
crack initiation strength {the tension required to initiate cracking): and the 
crack propagation strength (the tension below which the crack propagation 
halts). For the various eases described below, we specify the values adopted 
for each of these parameters. 


4 Effects of Obliquity on Cycloids 

4 1 Fitting Global Cycloid Pattern!: 

Figs. 8 and 9 show hypo! hetirul cycloids formed when !). i° of obliquity is added 
to the model using an initiation stress of 68 kPa. a propagation strength of 
38 kPa and an effective propagation speed of 4.8 ktn/hr (all typical param- 
eter choices that give generally rwsouaMe agreement with observed cycloid 
shapes). Figs. 8 and 9 assume westward and eastward propagation, respec- 
tively. As shown in the separate parts of Figs. 8 and !), the pattern depends 
on where in the orbit, relative to pertjove, Jupiter crosses Europa’s equator 
(or equivalently, the angle from the ascending node to the position of perksen- 
ter, known in celestial mcchanicH as argument of periccuter, m). For example, 
Fig. 8b shows the hypothetical cycloid pattern when Jupiter crosses Europa’s 
equator 1/8 of an orbit after perijove (i.e. when u —■ 3.15°). and so forth. 

An obliquity of 0. 1" is enough to break the equatorial .symmetry seen in Fig 1. 
For example, in Fig. 8a and Fig, 9a. when the argument of pcriccntcr is zero, 
“boxy” cycloids just west of 0“ longitude are shifted southward while ISO” 
away they are shifted northward. Moreover, some cycloids propagate, across 
the equator. 

Depending on the argument of periccuter. at. the shift in the hypothetical 
cycloid pattern changes. With st - 0“ (Fig. 8a). cycloid patterns near the 
equator around 30“ longitude are shifted south ~ 6°, but with ui — 315" (Fig. 
8b) cycloids in at this same region arc shifted north by ~ 4" in latitude 
With = 270” and 225° (l'ig. 8c & d). these cycloids are shifted even farther 



modeling the formation of cycloids can still uniquely constrain formation lon- 
gitudes As long as the location of cycloids have been moved in longitude only 
since their formation, preserving their formation orientations, fits are unique 
as statistically significant. 

5 Discussion 

We have shown the effects of obliquity on the formation of hypothetical cy- 
cloids on Europe. Obliquity introduces a component of stress due to the daily 
migration of the tidal bulge in latitude, breaking the symmetry in idealized 
crack patterns about the equator. These effects may account for many of the 
previously unexplained characteristics of the observed cycloid pattern on Eu- 
rope's surface- Cycloids that cross the equator arc not symmetrical, and the 
“boxy” cracks just west of 180“ arc south of the equator, while those just west 
of 0" arc north of the equator. 

Adding obliquity results in cycloids propagating across the equator without 
changing their eastward or westward sense of propagation, now explaining 
how such observed cycloids can form. The model now predicts long chains of 
equator-crossing cycloids near the same longitudes as the "boxy” patterns. Ft 
is, in fact, the north-south shifting of these regions that give rise to the cycloids 
that cross the equator. Observations of equator-crossing cycloids match such 
concentrations, near the Argadnel Regio region and ISO" in longitude away 

“Boxy" or angular cycloids resemble the observed crack pattern in the Ar~ 
gadncl Regio region, west of the subjovian point (near 210° longitude) (Bart 
Kt of., 2003). Without obliquity, models produce such shapes at this longitude. 


but they predict a symmetry about the equator. However, observed features 
are concentrated south of the equator and centered around a latitude of 5”S. 
Maybe the conditions needed to form angular cycloids north of the equator 
did not exist when cycloids in this region formed, naturally producing the 
observed asymmetry, but the addition of the stress from the obliquity into 
the model of cycloid formation lends to an asymmetry about the equator 
which can explain the observed asymmetry. A tenth of a degree of obliquity, 
the amount estimated by Bills {2(106), is adequate to shift the center of the 
stair-shaped cycloids south by at least 5° in the Argadnel Regio region, if these 
cracks formed when the argument of pericentcr, w was between 1S0° and .115“. 
Larger amounts of obliquity would produce larger shifts in latitude. 

Moreover, angular cycloids in the Argadnel Regio region have prominent seg- 
ments with orientations that are north-cast/south-west or north-west/south- 
east. However, cycloid formation models without, obliquity would produce an- 
gular cycloid with umth-west/snuth-east orientation south of the equator. The 
asymmetry introduced by the inclusion of stress from obliquity in the model 
of cycloid formation also allows angular cracks with north-east/south-west 
orientations to be produced south of the equator. Thus, more completely re- 
producing the angular cycloids observed in the Argadnel Regio region. 

As we have shown, the shift in the cycloid pattern can be to the north or 
the south depending on when in the orbit Jupiter crosses Europe's equator. 
Formation of the crack pattern in this Argadnel Regio region, being south of 
the equator, would require the argument of pcricenter to be between IH0* 
and 315“. Similarly, the “boxy" cycloids in the northern hemisphere near (T 
longitude would also require the argument of pericentcr to be between 180" 
and 315°. Moreover, equator-crossing cycloids between 0” and 30“ longitude 



specific cycloids may be able to place a more precise upper limit, to the size of 
Europa\s obliquity (Sarid el aL, 2006). 

Hoppa et al. ( 1999a) generated predictions of the sense of strike-slip motion on 
Europe's surface (i.e. left -lateral versus right-lateral), using the surface stress 
due to diurnal variations in stress produced by Europe’s eccentricity. In order 
to compare Hoppa s model to observations, Sarid ef at. (2002) produced maps 
of the observed strike-slip motion seen along Europan ridges, and attributed 
mismatch between the model and observations to polar wander. That work 
should mow be revisited, taking into account the effects of obliquity, which 
may be an alternative way to match the model anti observations. 

Recent additions to cycloid formation models focus on cycloid cusps (Green- 
leer and Katteuhorn, 2008). In these models normal and shear stresses are 
resolved onto the tip of the cycloidal crude while it remains dormant. Shear 
motion at a critical point in the orbit leads to the formation of a ridlerack 
which propagates away from the cusp to become the next, arc in the cycloid 
cha in. Cusp angles will be affected by the ratio of shear to normal stress needed 
to produce the next arcuate segment and will subtly affect characteristics of 
the cycloid. However, these additions can not explain the formation of equator 
crossing cycloids (Fig. 5). 

fi Conclusions 

Europa’s obliquity may haw played a substantia! tote in defining crack pal 
terns. The observed north-south asymmetry of cycloids may be evidence for 
the small but finite obliquity of Etiropa predicted by Bills (2005). Obliquity 


may be the key to producing aspects of the observed cycloid distribution such 
as equator crossing cycloids and an asymmetry in the "boxy” patterns. All at 
these improvements of the fit between observation and theory seem to require 
that the argument of perk-enter was between 180” and 345* at the t ime of crack 
formation, suggestive of complex dynamics, such as a linkage between the as 
raiding node and the orbital resonance. Future modeling of cracks should no! 
assume that the obliquity is zero or has no effect cm tectonic formation: While 
couch more work i.s needed to understand the role of obliquity in the tectonics, 
we have shown that, it was probably strong enough to fn; included in cycloid 
and stress models. 
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a Westward Propagating Cycloids: Sf NSR. 0" ObJiquit} 
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Fig. 1. Maps of hvpot hi’tioa! cycloids fcwiruxi using typical parameters found from 
modeling cycloids near the south pole: The initiation stress is 68 kPa, the propaga- 
tion strength is dS ki>a and the propagation velocity is 4.8 krn/fu {Hurford ft. a L 
2007). Only diurnal stress caused by Europa's eccentricity contributes to the stress 
held in which these hypothetical cycloids form, (a) shows westward propagating 
cycloids while (b) shows cycloids which propagate towards the oast- 
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Fig. 4, An example of a cycloid which trends mostly north sonth as it crosses the 
equator. Most equatorial crossing cycloids do not share this diamcknist ir. 
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Fig. 7. Obliquity causes # diurnal migration of tho tidal bulge in latitude. At some 
point during Europe's orbit. Jupiter U> aligned with the equatorial plain* of Europa 
{a) but obliquity causes Jupiter to move north of the equatorial plane until it madias 
ilsi farthest, northern point 1/4 an orbit later (b). Then Jupiter moves* south until 
it once again Is aligned with the? equatorial plane (c) iwifore moving to is farthest 
southern point 3/4 of an orbit later (d). 
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Fig. 8- Maps of cycloids that propagate westward, using typical parameters found 
from modeling cycloids near the south pole (Hurford d td. 2007). Only diurnal 
stresses caused by Europa a eccentricity and 0.1° of obliquity contribute to the stress 
field in which these hypothetical cycloids form. Each panel (a-h) shows 11 different 
perioentsr position measured relative to the ascending node. i.e. tho argument of 
pcrioentcr. in increments of 45”. 
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Fig. 8. [Continued) 
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Fig. 8. { Contiri ual } 
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Fig. 8, ( Continued) 


45 


Eiastwtird Propagating Cycloids 0.1° Obliquity, Argument of Pert'center 0“ 
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Fig. ff. Maps of cycloids that, propagate (vest ward, using typical parameters found 
from modeling cycloids near the south polo (Hurford tl til. 2007). Only diurnal 
stresses caused by Emopics eccentricity and 0.1' : of obliquity contribute to t lie stress 
hold m which these hypothetical cycloids form. Each panel (a h) shows a different 
fH'iieenter position measured relative to the ascending node. he. the argument of 
[toficiiiiter. in increments of 45®. 
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Fig. 9. ( Cim tinned ) 
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Fig. 0. (Continued} 
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Fit witii 70" longitude shift, 0,25° obliquity, to = 345 " 
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Fig. 10. A preliminary lit. of an equator-crossing eydokh which is shown in Fig, 
5. The cycloid is fit using I he method outlined in Harford et «/, (2007) and the 
values of the various parameters (clank speed, initiation strength and propagation 
strength) arc listed below the figure 
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